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Soft-collinear effective theory (SCET) is used to sum Sudakov double-logarithms in the a; — > 1 
endpoint region for the deep inelastic scattering structure function. The calculations are done in both 
the target rest frame and the Breit frame. The separation of scales in the effective theory implies 
that the anomalous dimension of the SCET current is linear in In /x, and the anomalous dimension 
for the A'^*'' moment of the structure function is linear in InN , to all orders in perturbation theory. 
The SCET formulation is shown to be free of Landau pole singularities. Some important differences 
between the deep inelastic structure function and the shape function in B decay are discussed. 
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I. INTRODUCTION 

The deep inelastic scattering cross-section is the in- 
clusive cross-section for lepton scattering off a hadronic 
target at large momentum transfer. The cross-section is 
conventionally written in terms of structure functions of 
the momentum transfer and a dimensionless variable 
< a; < 1. The structure function(s) F{x, Q'^) cannot be 
computed in perturbation theory, but its dependence 
can. The structure function contains large logarithms 
of the form {cts IhQ^/Aq^q) , which can be summed by 
evolving F{x, Q^) from the large scale to a lower scale 
/i using the renormalization group equations, /i is chosen 
to be a few GeV, paramctrically of the order of Aqcd, 
but still large enough that perturbation theory is valid. 

As a; ^ 1, there are additional large logarithms that 
need to be summed to get a reliable evaluation of the 
scattering cross-section. The invariant mass of the final 
hadronic state is 
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as X 



1. The total cross-section is in- 
frared finite, even though the real and virtual emission 
processes are separately infrared divergent. In the re- 
gion X ^ 1, real gluon emission is suppressed, and the 
cancellation between real and virtual emission becomes 
more delicate, leading to large corrections to the cross- 
section. The form of the perturbation series is most con- 
veniently described in moment space, where x — > 1 cor- 
responds to large moments, N — > oo, with the heuristic 
rule 1 — a;~l/A^. AsA^^oo, the structure func- 
tion moments contain terms of the form a!! In^ N with 
s < 2r. These Sudakov double-logarithms are important 
in the endpoint region. The summation of these terms 
is well-known, and has been discussed extensively in the 
literature Q . The general result is that the A^'^ moment 
of the structure function at Q^, Fm{Q^) can be written 

as [an 

lnF^(Q2) ^ fo{as\nN)\nN + h{as\nN) 

+asf2{aslnN) + ... . (2) 

The exponential of /o («« In A^ ) In A^ gives the leading Su- 
dakov double-logarithmic series. 



In this paper, Sudakov double logarithms in the end- 
point region are calculated using soft-collinear effective 
theory (SCET) HH- SCET allows one to compute the 
cross-section in the endpoint region in a systematic ex- 
pansion to any desired accuracy. The results are free 
of Landau pole singularities. The calculations are de- 
scribed in detail in both the target rest frame and the 
Breit frame, and give an instructive example of the use 
of SCET. There are several unusual aspects of SCET 
which are discussed here; the frame dependence of the 
way in which infrared divergences cancel between the soft 
and coUinear modes, and the relation between structure 
functions and local operators which is different for deep 
inelastic scattering and B decays. Consistency of the ef- 
fective theory implies that the SCET anomalous dimen- 
sion is linear in In jjb, which leads to the form Eq. (0) of 
the perturbation series. The SCET calculation to the ac- 
curacy presented here gives the moments of the structure 
function including the first two exponentiated series /o,i 
in Eq. (O, as well as all terms of order as which do not 
vanish as A^ ^ oo. The method used parallels that for 
B in Ref. 0. Deep inelastic scattering structure 

functions in the Breit frame (but not in the endpoint re- 
gion) were considered in Ref. 0. The SCET formalism 
used in this paper is described in Refs. Q. 



II. OUTLINE OF CALCULATION 

The calculation of the deep inelastic scattering cross- 
section will be performed using a sequence of effective 
field theories. The scattering amplitude involves the in- 
teraction of a lepton beam with a hadron target via a 
virtual photon. The leptonic interactions are calculable 
using QED, and will not be discussed here. The quantity 
of interest is the interaction of the virtual photon with 
the hadronic target. 

At scales much larger than Q^, the interaction of pho- 
tons with hadrons is described using the electromag- 
netic current of quarks interacting via the full QCD La- 
grangian. The hadronic scattering amplitude is the ma- 
trix element of the electromagnetic current between the 
initial and final hadronic states (see Fig. 
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FIG. 1: The deep inelastic scattering process. The incom- 
ing lepton with energy E scatters off a hadronic target with 
momentum p to produce the final hadronic state X. 

= {X\f\P) . (3) 

The cross-section involves the square of the scattering 
amplitude, and can be written in terms of the product of 
two currents, summed over intermediate states, 

J2{P\j^\X){X\f\P) . (4) 

X 

The invariant mass of the final hadronic state for generic 
values of x is of order Q^. The final hadronic state X 
can be integrated out at the scale Q^, and the product 
of currents in Eq. Q is replaced by a sum over local 
twist- two operators when one takes moments of the cross- 
section. This is the conventional method of computing 
the deep inelastic scattering cross-section. 

As X — + 1, the invariant mass of the final hadronic state 
tends to zero. The invariant mass of the final hadronic 
state is taken to be of order Q^X, with ^ Q^X ^ 
-^QCDi ^'^^ results are computed in an expansion in 
A. More details about this power counting scheme are 
given in Sec. lIII CI The scale Q^X is an infrared scale for 
the theory at Q, and can be set to zero at leading order in 
an expansion in A. In this limit, the final hadronic state is 
massless, and cannot be integrated out. Instead, the final 
hadronic state can be treated as a massless light-like jet, 
and is described by a coUinear quark in SCET. Thus for 
X ^ 1, the QCD current at scale Q is matched onto an 
SCET current at scale Q. The coordinate axes are chosen 
so that the outgoing quark travels in the = (1, 0, 0, 1) 
direction, and is described by the coUinear field ^„ of 
SCET. In the target rest frame, the incoming quark has 
momentum of order the hadron target momentum Aqcd , 
and is an ultrasoft quark. In the Breit frame, the struck 
quark is back-scattered, so the incoming quark is a n 
coUinear quark, and is described by the field S,n, where 
fi = (1,0,0,-1). The first step of the calculation is the 
matching of the electromagnetic current ip^f^tp onto the 
SCET current ^nl^'^ or ^n'y'^S.n in the target or Breit 
frame, respectively. The matching coefficient is evaluated 
in Sec. |IV| 

The next scale in the problem is the invariant mass 
of the hadronic state, pj^ ~ Q^A ^ Q'^{1 — x). The 



SCET current is run from /i = Q to this scale using the 
anomalous dimension, which is computed in Sec. [V] in 
both the target and Breit frames. 

At the scale (5^(1 — x) ^ Q^X, p\ treated as large, 
and so the final state can be integrated out. The time- 
ordered product of two SCET currents can be replaced 
by a bilocal light-cone operator whose target matrix ele- 
ment is the parton distribution function. This is done by 
integrating out the ^„ field in SCET; the resulting oper- 
ator is written in terms of ultrasoft quark fields ipu in the 
target frame, and in terms of fi-collinear quark fields in 
the Breit frame. The bilocal operator is closely related 
to the CoUins-Soper operator 0]. The matching at scale 
Q^A is discussed in Sec. IVII 

The last step, discussed in Sec. lVIII is to run the bilocal 
operators from Q^{1 — x) to some low scale /i, and match 
onto local twist-two operators which give the moment 
sum rules for the deep inelastic structure functions. 

The computations in this paper are given in Feynman 
gauge. The results, however, are gauge invariant, and 
valid in any gauge. In the effective theory, one has sep- 
arate gauge invariance for the ultrasoft, n-coUinear and 
n-collinear gluons. This has been checked by explicit 
computation. 

The entire analysis is presented for QCD with a sin- 
gle quark flavor of unit charge, to avoid unnecessary in- 
dices. The final answer is given by summing the results 
of this paper over all flavors weighted with the square of 
their electromagnetic charges. In addition to the currents 
^nl^i^ and (,nl^S,n, one also has the hermitian conjugate 
currents "07^^" and ^nl^£,m which give the antiquark 
contribution, or equivalently, the crossed-graph contri- 
butions. By charge conjugation invariance, the matching 
coefficients and anomalous dimensions are the same for 
quarks and antiquarks. The final result thus has a sum 
over both quark and antiquark distributions. 

III. KINEMATICS 

The scattering process is e" +p +X. The proton 
momentum is P, the incoming momentum of the virtual 
photon is q, the momentum transfer is = —q^ 3> 0, 
and X is defined by 



The coordinate axes are chosen so that the virtual photon 
is in the z direction. It is useful to introduce the null 
vectors n'' — (1, 0, 0, 1) and — (1, 0, 0, —1) which point 
in the z and — z direction, respectively. They satisfy the 
relations ~ 0, = 0, fi ■ n = 2. Any four-vector 
can be written as 

a" = ia+n^-^ia-n^ + a'^, (6) 

where 

= n ■ a , a~ = n ■ a , (7) 
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and a_L is in the x — y plane. The dot product of two 
four-vectors is 

a-b ^ ]j^a+b^ + ^a^^"'' - • , (8) 

and the integration measure can be written in light-cone 
coordinates as 

d'^fc = i d/c+ dfc- d'^~2fc_L . (9) 
A. Target Rest Frame 



so that 

1-. = ^, (14) 

and 

p\ = Ql^- (15) 

The ± components of momentum in the Breit frame 
are given by multiplying the ± components of momentum 
in the target rest frame by {Q/xP~^)^ ~ (Q/P^)^ for 
a; « 1. 



In the target rest frame, qj^ = 0, and = —q^ = 

—q^q^. The deep inelastic limit — > oo with x fixed is 

given by taking g~ ^ oo at fixed so that q^ ^ q^ . 
Then 



where P+ = Mt is the target mass, so that 













xP+ ' 


v\ 


= P+(1-.t) , 


Vx 


= p-+q-^q 


v\ 





(11) 



where px = P + g is the momentum of the final hadronic 
state. 



B. Breit Frame 

The virtual photon carries only momentum, and no 
energy in the Breit frame. The Breit frame is obtained 
from the target rest frame by boosting along the 2;-axis, 
so that the proton and virtual photon have no _L com- 
ponent of momentum in either frame. The momentum 
components in the Breit frame are 



g+ = -Q, 
(f = Q, 
P+ = Q + V 

= 

Px = 

Px = Q + i' 



where are fixed by setting P^ fa Ql = Mt and 

Q ^ Q 

^ " Q + 1+- MtIQ ^ Q + 1+' 



(12) 



(13) 



C. Power Counting 

The endpoint region is ^ Q^{1 ~ x) ~ Q^A, with 
^ Q^X ^ ^QCD- this region, the final state still 
involves a sum over many hadronic states, but has small 
invariant mass and is jet-like. The dimensionless power 
counting parameter A ^ 1 is introduced as the expansion 
parameter, and 1 — a; ^ A. 

The power counting is simplest in the Breit frame, 
where ~ QX and £~ ^ Aqcdj so that P"*" ~ Q, 
P- ^ Aqcd, Px ~ <3A, Px ^ Q, q+ ^ Q, q^ ^ Q. 
There are three important scales: (1) Q^, the invariant 
mass of the virtual photon, (2) Q^X, the invariant mass of 
the final state hadronic jet, and (3) ^ ^qcd^ the in- 
variant mass of the target. The scale Q^X^ does not play 
an important role in deep inelastic scattering; it does for 
the shape function in B decays 0- 

Particles with p^ ^ Q, p'^ ^ QX^ and pj^ ^ QX travel 
in the n direction, and are described by 7i-collinear fields 
^np- Pi(^) SCET. The large components of momen- 
tum, and are explicit labels on the field, and mo- 
mentum of order QX^ is the Fourier transform of the co- 
ordinate X. This is analogous to use of label-momentum 
for non-relativistic quarks in NRQCD '3| . Similarly, par- 
ticles with momenta p'^ ^ Q, p^ ^ QX^ and pj_ ~ QX 
travel in the n direction, and are described by n-collinear 
fields p^. Particles with momenta of order QX^ are 
described by ultrasoft fields. In the Breit frame, the out- 
going quark is described by a n-collinear field, and the 
incoming quark is described by a n-collinear field. The 
choice of coordinate axes is such that the _L components 
of label momentum are zero. The fields will be referred 
to as for simplicity. 

The Breit frame is the natural frame to use to describe 
deep inelastic scattering near x — \. The power counting 
automatically implies that 1 — x — > 0, by Eq. (|14|l . Nev- 
ertheless, it is instructive to also give results in the target 
rest frame. The target frame is the best frame to com- 
pare deep inelastic scattering with B Xs^- The target 
rest frame is also the natural frame to use for generic 
values of x, and is a ^-independent frame. The boost 
to the Breit frame depends on a;, though for a; « 1, the 
boost factor is approximately constant. In the target 
rest frame, the incoming quark has momentum of order 
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FIG. 2: One-loop vertex correction to the electromagnetic 
current in QCD. 



Aqcd, and is described by an ultrasoft field. The outgo- 
ing quark has momentum components pj^ ~ AAqcd/Q, 
~ O^/^QCD, ~ QA. The outgoing quark turns 
into a jet moving in the n direction, and so the outgoing 
particle can be described by a n-collinear field. 



IV. MATCHING THE CURRENT AT 
QCD TO SCET 



FROM 



The electromagnetic current in QCD is matched onto 
the SCET current at the scale Q. The QCD current is the 
operator and the SCET current is fnW„7''W^^ft 

in the Breit frame, and fnW^re7^'*/' iii the target rest frame. 
Here Wn,n are coUinear Wilson lines which are required 
by coUinear gauge invariance. The matching condition 
will be computed in (a) pure dimensional regularization, 
i.e. using dimensional regularization to regulate both the 
ultraviolet and infrared divergences, and (b) by using di- 
mensional regularization for the ultraviolet divergences 
and ofF-shellness for the infrared divergences. 

The one-loop vertex graph for the electromagnetic cur- 
rent in QCD is shown in Fig. [3 where pi is the incoming 
quark momentum, and p2 = pi -I- g is the outgoing quark 
momentum. The QCD one-loop graph gives 



1 



(16) 



where Cf = 4/3 is the Casimir of the fundamental rep- 
resentation. 

We first consider the computation in pure dimensional 
regularization, which greatly simplifies the computation 
of matching conditions in effective field theories. In pure 
dimensional regularization, the matching coefficient is 
obtained by computing the finite parts of on-shell di- 
agrams and dropping all the 1/e terms, regardless of 
whether they arise from ultraviolet or infrared diver- 
gences y\. The reason this procedure works is that the 
ultraviolet divergences in the full and effective theories 
are canceled by the counterterms in the respective theo- 
ries. The remaining 1/e terms are infrared divergences, 
which must agree between the full and effective theory. 
The 1/e terms cancel in the matching condition, which is 
the difference between the full and effective theory. Thus 
the matching condition is the difference of the finite parts 
of the full and effective theory computation. There is one 



additional simplification — on-shell graphs in the effec- 
tive theory are usually scaleless integrals, which vanish 
in pure dimensional regularization, and so have no finite 
part This eliminates the need to compute the effec- 
tive theory graphs to determine the matching condition, 
which is given by the finite part of the full theory graphs. 

We will compute Eq. H16|) keeping the 1/e terms to 
compare with the matching computation using an in- 
frared regulator. The incoming and outgoing quarks in 
Fig. 121 have invariant masses that vanish in the limit 
A 0, so the matching coefficient is obtained by evalu- 
ating the graph on-shell with = p2 _ Evaluating 
the integral in d = 4 — 2e dimensions gives 
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21n^ 



eiR 



? ? 9 



(17) 



where we have distinguished the infrared and ultraviolet 
divergences by the subscript on e. However, it is impor- 
tant to keep in mind that all e's are equal. The integral 
has a 1/ejp, infrared divergence arising from a combina- 
tion of soft and coUinear divergences. It is this double 
divergence that leads to the Sudakov double-logarithmic 
behavior in the endpoint region. In pure dimensional 
regularization, the wavefunction graphs are scaleless. 



47r 



1 



1 

eiR 



(18) 



and vanish. The net on-shell matrix element of the elec- 
tromagnetic current in the full theory is the difference of 
Eqs. H17I18|1 plus the counterterms, which gives (includ- 
ing the tree-graph) 



{P2\j'" \P1 



1 + 

I 4r 



21n^ 



eiR 



1 2 A*^ Ql fi 




hc.t. . 

(19) 

The 1/euv terms cancel, so there is no counterterm. This 
cancellation is required, since the electromagnetic current 
is a conserved current and has no anomalous dimension in 
QCD. The graphs in the effective theory are all scaleless, 
and vanish in dimensional regularization. The matching 
coefficient of the current in the effective theory is the 
finite part of Eq. ((T^ . 



Cifi) = 1 



47r 



9 9 

1 2 A* q , M 



(20) 



The 1 /eiR terms in Eq. (|19|l . which are the negative of the 
1/euv terms in the effective theory, give the anomalous 
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dimension of the current in the effective theory, as we 
will see in the next section. 

The logarithms in the matching coefficient C(/i) can 
be minimized by choosing the matching scale ^ — Q, ed 
which 



c{Q) = i + ^fMcp 



An 



(21) 



The matching computation can be repeated by reg- 
ulating the infrared divergence by using off-shell initial 
and final states, with Pi = 7^ 0- The graph in Fig. [21 
gives 



V 



in 



1^ ~ " 2 In — In — 

euv M Q Q 



-'4 -"4 -'4 



(22) 



where the 1/e term is purely an ultraviolet divergence, 
since the infrared divergences have been regulated by the 
off-shellness. The evaluation of Eq. H22(l is considerably 
more complicated than that of Eq. H17|l . The wavefunc- 
tion graph is 



1 

euv 



1 - In ■ 



(23) 



so that the matrix element in the full theory including 
the tree-graph is 



(P2| \P1 




1 1 P2 



(24) 



The ultraviolet counterterm vanishes as before, as it 
must, since it does not depend on the choice of infrared 
regulator. The matching condition is given by subtract- 
ing from Eq. H24I) the matrix element in the effective the- 
ory. The effective theory integrals are no longer scaleless, 
since they depend onpf, and must be evaluated to obtain 
the matching condition if an off-shellness is used to regu- 
lated the infrared divergence. The matrix element in the 
effective theory with an off-shellness is given in the next 
section, where the anomalous dimension of the SCET 
current is computed. Taking the difference of the result, 
Eq. (|34|) . and Eq. H24|) gives the same matching condi- 
tion as before, Eq. (^0)) . The computation of Eq. is 
clearly simpler using dimensional regularization to regu- 
late the infrared divergence, since it does not require the 
effective theory computation at this stage, leaving it to 
the next section where it properly belongs. 



V. ANOMALOUS DIMENSION OF THE SCET 
CURRENT 

The electromagnetic current in QCD is matched at the 
scale Q onto the SCET current. In the Breit frame, the 
current is 

.r = CiQ) ^n,QWnrWkn,Q, (25) 

where the n-collinear quark has label momentum n ■ p ^ 
Q, p± — 0, and the fi-collinear quark has label momen- 
tum n ■ p = Q, p± = 0. C{Q) is the matching coefficient 
computed in Eq. |2U. 

In the target rest frame, the current is 

,f = CiQ) L.,-Wnr^u, (26) 

where ^„ is an ultrasoft quark, Wn is a collinear Wilson 
line, and ^n,q- is a n-collinear quark field with labels 
n ■ p = , p± = 0. 



A. Breit Frame 

The one-loop anomalous dimension of the SCET cur- 
rent in the Breit frame is given by the graphs in Figs.Ol 
as well as wavefunction renormalization graphs. 

The ultrasoft graph, Fig.l^fa), gives 
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(27) 



Doing the fc+ integral by contours and using the substi- 
tution — xp^ gives 



dx 



2n J {2n) 



d-2 



8n 



[1-X + iO+] [ptPiX -k]^+ iO+] 
-^^r(6) [-p+P^x]- 
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The n-collinear gluon graph Fig. Ol^b) gives 

2^ 2e / d'^k fin" fin - {p2 - k) 



In = -ig Cpf^ 
1 

X ■ 



{2n)'^ 2 2{p2-kf+i0+ 
1 



(28) 



r 



fl-k + iO+ fc2 + iO+ 

d'^k n ■ {p2 — k) 



-2ig''CFfJ. 
1 



{2ny {p2 - fc)2 + iO+ 
1 



(29) 
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FIG. 3: One loop correction to the electromagnetic vertex in the Breit frame from (a) ultrasoft, (b) n-coUinear and (c) n-coUinear 
gluons. 



Evaluating the integral by contours, and letting k — 
gives 



1 dz 
2^ 
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(30) 



The rt-collinear gluon graph Fig. |3[c) gives 



fin - {pi 
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Tl 

(2Try -n-k + iO+ 
k) flfia 1 



' 2(pi - fc)2 + iO+ 2 fc2 + iO+ 



(31) 



which is Eq. ipn|l with p| 
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The remaining graphs are the wavefunction graphs. 
The ultrasoft gluon contribution to wavefunction renor- 
malization vanishes in Feynman gauge, since n2 = n2 = 
0. The collinear wavefunction renormalization graph is 
the same as the massless quark wavefunction renormal- 
ization in QCD, Eq. H23|l . since the interaction of n- 
collinear quarks with rt-coUincar gluons is the same as 
the interaction of quarks with gluons in full QCD. 

The matrix element of the current in the effective the- 
ory is given by the sum of Eqs. H28I30I32|I and subtracting 
half the wavefunction renormalization Eq. H23() for each 
external quark. The net result is 
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(33) 



The ultraviolet divergence of this result agree with the ul- 
traviolet divergence in the effective theory inferred from 
Eq. p9|) . Exactly on-shell, the effective theory integrals 
are scaleless and vanish because the 1 / euv ultraviolet di- 
vergences cancel the 1/eiR infrared divergences. The ul- 
traviolet divergence in the effective theory should there- 
fore be the negative of the 1 / eiR terms in Eq. (|19|l , which 
agrees with the divergence in Eq. The SCET cur- 

rent operator is multiplicatively renormalized, and there 
is no operator mixing. This allows one to compute the 
anomalous dimension of the SCET current directly from 
the 1/eiR terms in the full theory matrix element. In 
cases with operator mixing, the 1/eiR terms in the full 
theory matrix element gives the value of the anomalous 
dimension matrix times the operators coefficients evalu- 
ated at the matching scale. 

The infinite part of the matrix element is canceled by 
the vertex and wavefunction counterterms in the effective 
theory so that the renormalized matrix element is 



{P2\f\pi),, 



47r 



-2 In- 



2 In- 



9 

-Pi 



1 

m 



2 In 



-pi 



In- 



9 

-P2 



In- 



'P\ 



In- 



9 

-Pj. 

^2 



-pI, 

_ 5^ 
6 

(34) 



The infrared divergence as ^ in the full theory cal- 
culation, Eq. (|24|) agrees with the infrared divergence of 
the effective theory calculation, Eq. (I34II . and the differ- 
ence gives the matching condition, Eq. (|20|1 . which is free 
of infrared divergences, and depends only on Q^, not on 

pi- 

The ultraviolet counterterm for the ultrasoft graph 
Fig. OJa), ri-coUinear graph Fig. Olb) and n-collinear 
graph Fig. OIc) and wavefunction graph are from 

Eqs. (I:^8l3()l3:^l:^3li 



ultrasoft: -r-Cp 
An 



n-collinear: — Cp 
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n-collinear: 
wavefunction: 
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respectively. The individual counterterms are sensitive 
to the infrared through their dependence on the small 
components of momentum, and p J . The total coun- 
terterm is the sum of the four terms above, 
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and depends only on the label momenta "p^ and p\ which 
are both Q. The counterterms Eq. H36|l give the anoma- 
lous dimension for the coefficient of the current in the 
effective theory. 
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d/^ 
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The SCET current anomalous dimension depends on In /i, 
because the one- loop diagrams have terms from 

combined coUinear and soft divergences. Consistency of 
the effective theory implies that to all orders, the anoma- 
lous dimension is at most linear in In /i, as will be shown 
in Sec. |IX| 

The cancellation of 1/e Inpj^ and 1/e Inp^ between the 
soft and coUinear graphs might suggest that the coef- 
ficients in the two sectors must be the same, i.e. that 
the two contributions must have the same value of as- 
This suggests that SCET should use two-stage running, 
in which all coupling constants are evaluated at a single 
scale /i, unlike NRQCD, which requires one-stage run- 
ning using the velocity renormalization group [8LllO|. In 
NRQCD, there is a cancellation of infrared divergences 
between the soft and ultrasoft sectors that naively sug- 
gests that both should have the same value of as- How- 
ever, this is not the case, and a proper treatment of 
NRQCD has the soft coupling constant evaluated at the 
scale mv and ultrasoft coupling constant evaluated at 
the scale mv^ ^| . It has been pointed out that one- 
stage and two-stage running give the same result in SCET 
for quantities which have been computed so far llj. In 
NRQCD, the difference between one- and two-stage run- 
ning first occurs at order in the power counting. SCET 
anomalous dimensions have so far been computed only to 
leading order in A, and do not distinguish between one- 
stage and two-stage running. 




P2 



Pi 



(b) 



FIG. 4: One loop correction to the electromagnetic vertex in 
the target rest frame from (a) ultrasoft and (b) n-coUinear 
gluons. 



B. Target Rest Frame 

In the target rest frame, the electromagnetic current in 
the effective theory contains a n-collinear quark and an 
ultrasoft quark. The diagrams in the effective theory in 
the target rest frame are those in Fig.Qland wavefunction 
graphs. Note that there is no n-collinear graph, since the 
current Eq. H25|l does not contain Wn in the target rest 
frame. 

The n-coUinear graph Fig. 2{b) is the same as in the 
Breit frame, and is given by Eq. (|30|l . The result is frame- 
independent, since the final expression only depends on 
P21 not on the individual components or the momen- 
tum -px of the incoming quark. The non-zero wavefunc- 
tion renormalization graphs also agree between the two 
theories. The n-collinear wavefunction renormalization is 
the same as the corresponding graph in the Breit frame, 
and the ultrasoft wavefunction renormalization of the ul- 
trasoft quark is the same as the n-collinear wavefunc- 
tion renormalization of the n quark, since both are equal 
to wavefunction renormalization of a massless quark in 
QCD. Each wavefunction graphs depends only on p^ for 
a single particle. 

The remaining graph is the ultrasoft graph Fig. Qfa) . 
The ultrasoft graph in the target rest frame must contain 
the contributions of both the ultrasoft and fi-collinear 
graphs in the Breit frame. The ultrasoft graph in the 
target rest frame is 



d'^fc 



1 



{2nY n- {p2- k)+ iO+ 

h-^ 1 



(pi - + iO+ ^" fc2 iO+ 

ig'CF^" 



{2n)'^ [n ■ {p2 - k) + iQ+]{pi ~ kf + iQ+ k^ + iO+ ' 

(38) 

Doing the fc+ integral, and using fc"*" = zp^ gives 
dz f d'^-^ki 



I 



7pt-ptz){z{l-z)pl-^_) 



8 



ipt-Pi^) 



(39) 



The ratio /p'^ is 



= 1 — X 

Pi 





2 2 2 

P2 ^ ^ ^ ^ Q(^x) . (40) target rest frame 



FIG. 5: Tree graph for the product of two currents in the 



ptP2 



In SCET, terms of order 1 — x are of order the expansion 
parameter A in the power counting. Equation H39|) can 
be evaluated in the hmit /p^ — s- 0, which simpUfies 
the computation, and gives 




This result is identical to the sum of the ultrasoft and n- 
coUinear graphs in the Breit frame, given in Eqs. (|28I32|) .-^ 
The ultrasoft graph in the target rest frame has no 
divergence. Since it is the sum of the ultrasoft and n- 
collinear graphs in the Breit frame, these graphs must 
have terms of opposite sign, which agrees with the 
explicit one- loop computation in Eqs. (|28I32|I . This can- 
cellation is expected to persist at higher orders. 

The 1/e term in Eq. H41(l depends on the infrared reg- 
ulator through using Eq. (|40|) . This infrared depen- 
dence is canceled by the n-coUinear graph, which can only 
depend on p| . The analog of Eq. (|35|l is 



VI. MATCHING AT Q^{1 ~ x) ONTO THE 
PARTON DISTRIBUTION FUNCTION 

At the scale Q^{1 — x), the invariant mass of the final 
hadronic state p^ can be treated as large, and the final 
state can be integrated out. This is done by integrat- 
ing out the n-collinear modes from the effective theory. 
Since the current Eq. H25I26|1 contains ^„ fields which are 
integrated out, one matches the product of two currents, 
rather than a single current, onto the effective theory be- 
low Q^{1 — x). The matching coefficients are determined 
by computing the matrix elements of 



271 



d4a;e''?-^j^(x)j''(0). 



(44) 



at fixed X and q^. Note that we have a product of cur- 
rents, rather than a time-ordered product. The matrix el- 
ement of the product of currents is given by the disconti- 
nuity in the matrix element of the time-ordered product.^ 
We will use this procedure, since time-ordered products 
can be computed using conventional Feynman diagram 
perturbation theory. 



ultrasoft: 
n-collinear: 
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so that the total counterterm is 



c.t. 
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(42) 



e pipi 



e ptP2 



(43) 



which is the same as Eq. ()36|l . and leads to the same 
anomalous dimension, Eq. (|37|l . 

The on-shell matrix element in the target rest frame 
has the same value as in the Breit frame, and leads to 
the same matching condition Eq. (|20|l . 



A. Target Rest Frame 

The tree graph for the product of two currents with 
integrated out is shown in Fig. |S1 The spin averaged 
matrix element of the tree level graph is 

27r 2 ^' 2 {p + q)^+ iQ+ ' 

= Disc ±T^^lliP-pLl 
271 [p + qY + iO+ 



T'"'5(l + g+/p+) 



where 



p^'n'' + p''n^' 
n ■ p 



Since p is at rest, oc = (1, 0, 0, 0) 



(45) 
(46) 

(47) 



^ The components of pi have different values in the two frames, ^ The discontinuity of a diagram containing iO"*" terms is defined 
but and pj/pf are equal, since the Breit frame is obtained by taking the difference of the diagram, and the diagram v^ith 

from the target frame by a boost along the z axis. iO"*" — > iO~ . 
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FIG. 6: Quark matrix element of the quark distribution op- 
erator. 



independent of the momentum of the target. 
Define the quark distribution operator by [ 



d^e-^^''^ 



(48) 



where Y is an eikonal Wilson line from to containing 
ultrasoft gluons A„, 



Y{n^,0) = Pexp 



n ■ Au{nz) dz 



(49) 



and 4'u are ultrasoft quark fields. The Feynman rules are 
given by taking the discontinuity of the diagram, since 
the operator is a product, not a time-ordered product. 
The quark distribution in a target T with momentum P 
is defined by 



f,/T{x) = {T,P\0,{xP+)\T,P) 



(50) 



The only difference between the quark distribution oper- 
ator Eq. (|48|l and the conventional CoUins-Soper defini- 
tion is the replacement of the full theory quark field by 
the ultrasoft quark field in SCET. 

The spin-averaged tree-level matrix element of the 
quark distribution operator is (see Fig. Ej) 

I i 1 , , 
Disc : --Tt i4 



Disc 



47r n • (p — r) 2 
I 2ip^ 



(51) 



(52) 



47r p+ - r+ + iO+ 
= <5(l-r+/p+), 

so the tree level relation is 

W^'"' = T'"'Og(-g+) . 

This is an operator relation independent of the matrix 
element since T'"' in Eq. (|47|l can be written in a form 
independent of the target state. The minus sign in the 
argument of Oq arises because q is an incoming momen- 
tum and r is an outgoing momentum. Equation (|52|l can 
also be written as a convolution 

W^" = T^"" I —S { 1 ^ 1 Og{r+) . (53) 



ft is convenient to write 
that Eq. 153(1 becomes 



— ^yP^ J and r+ ~ wp'^ so 



-yp 



dw . 



\ wJ 



(54) 



Note that y and w in this equation are defined with re- 
spect to the parton momentum p rather than the hadron 
momentum P. 

The analysis has been restricted to spin-independent 
structure functions for simplicity. It is straightforward 
to generalize this to spin-dependent structure functions, 
which involve spin-dependent quark distribution opera- 
tors HI. 

The one-loop graphs for the matrix element of the cur- 
rent product are shown in Figs. [7| Graphs [7Ja-d) have 
the same value in the theories above and below (5^(1 — x), 
since the interaction of an ultrasoft gluon with a coUinear 
quark is the same as the vertex generated by the Wil- 
son line Y in the operator Oq , and the external ultrasoft 
quark fields are unchanged at the matching scale. The 
matching correction is given byCfe-g); these graphs are 
absent below Q^{\ — x) since the n-collinear modes have 
been integrated out. 

The coUinear graph Fig. [Tfe) gives for the spin- 
averaged matrix element 

/el = Disc ^ / ;^TrV^";^^ + f 
2^7 {2ttY 2^' 2(p + g)2 

fin°- fifi ■ {p + q — k) ^ 1 _ 1 



2 2(j5 + g-fc)2 ' _f^.fc'^"/!c2 

2 + f d'^k n ■ {p + q) 
p^ 



Disc ^T''''n+ 



n- {p + q — k) 1 1 

{p + q — k)"^ —n ■ k 

Comparing with Eq. H29|l gives 



(55) 



Disc — 77^ p+T^ , 

27r V 87r2 y ^ (p + q)^ 



1 11, 

h -ln■ 
e^ e e 

\2 



Disc — ( 1 T'"' 
27r I 87r2 



^2 

2 



-{p + q? 

/i2 



1 



1 1 1 



In- 



(l-y) + zO+ 



1^^^2 Q^(?/~1)~»0+ 



VP- 



In- 



-2 + T2 



(56) 



where y > since g+ < 0, and = —Q^/q'^ > 0. The 
kinematic region where q~ < is infinitely far away in the 
effective theory, and is described by the effective theory 
with a different value of the label momentum. The graph 
with gluon attached to the other vertex gives the same 
contribution, Eq. H5t)|) . 
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The collinear graph Fig. [7ff) is given by the tree di- 
agram, times the negative of the wavefunction diagram 
Eq. evaluated with — > {p + qf , and gives 



Disc 



27r 



1 11, 



{p + g)2 + iO+ 
{p + q?^ 



1 



Disc — ) T^"^- 

27r V 87r2 / (1 - y) + iQ+ 



1 1 I ^^^ Q\y - I) - iQ+ 
2e 2 2 



0(y).(57) 



The graph Fig.^g) vanishes, since the coUinear gluon 
emission vertex is proportional to n", and = 0. 

The total collinear contribution is the sum of twice 
Eqs. (ESJ and Eq. (|57jl . 

^ Disc — ( ^\ T'^"^ \ 



^0^ 



In^^ 

7 ^2 
2+T 



3 ^ 2^^^ Q2(^_^) 
2e e ?//i2 



H2/-l)-zO+ , i^^^Q\y-l)-iQ+ 



yp 



(58) 



The collinear counterterms cancel the 1/e terms, so 
that the remaining matching contribution is finite. 
The collinear counterterm contribution (In— p2)/g = 
[ln(3^(y — l)]/f is no longer infrared sensitive, since the 
scale Q^{1 — y) is now considered a large scale. 



The discontinuity of the remaining terms can be writ- 
ten in terms of -|- distributions. The distribution 1/(1 — 
y)+ is defined by 



so that 



(1-2/)+ Jo 



1 



(1-2/)+ 



The discontinuity of 

i hi{y — 1 — irj) 



'a./lfifi, (59) 

(1-2/)+ 



(60) 



27r I — y + irj 



0{y) 



(61) 



is given by the difference of the expression for 77 0+ 
and 77 ^ 0^, 



„ i \n{y — 1 — im) , , 1 

^i^'^ — r — --^oiy) = 1 — 

ZTT 1 — y + IT] 1 — 2/ 



0(0 < y < 1). (62) 



The singular terms at y — 1 can be obtained by integrat- 
ing Eq. ifS^ over < y < A, where A > 1. Then 



2-K 



dy 



In (y — 1 — irj) 
1- y + irj 



in 



(63) 



irrespective of the sign of rj. This has no discontinuous 
part, so the integral of the discontinuity is zero. This 
gives 

Disc f i^fc^^0(y) = (64) 
27r 1-y + iT] (1-2/)+ ^ ' 
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Similarly, 



-J * 



In^jy - 1 - iij) 
1 — y + 17] 



-— [in^A- (-i7r)3sign77] , 
(65) 





so that 



Disc — 
27r 



dy 



ln^(y- 1 - ii]) 
I — y + irj 



(66) 



and 



Disc 



i ln^{y — 1 — ir]) 
27r 1 — y + ir] 



0{y) = -^'5(1-2/) 



21n(l-y) 



FIG. 8: Tree graph for the product of currents in the Breit 
frame. 



We have dropped terms of order 1 — y times the terms 
retained in Eq. (|68|) . All such terms have moments which 
vanish as iV ^ oo, so Eq. (|71|l is valid up to corrections 
which vanish as — > oo. Letting 



N = A^e'' 



1-y 



(72) 



(67) 



Using the above results gives for Eq. (|^ . 
Disc = Tf'Miy), 

"In(l-y)" 



finds 



M]s{M) 



2lT 



In" 3^1, 
^ Nfi^ 2 ^ Nfi^ 



TT 

T 



(73) 



M{y) 



^eio<y<iH2 



2 In 



■2 2 

-)2 



The logarithms are minimized if fj, = Q'^/N, so one can 
match from currents onto the bilocal operator at this 
scale using 



In 



2 3 g 



(i-y)H 



In- 



7 7r2 



-^(1 - y) 



MNiM) 



2^? 



TT 

y 



(74) 



(68) 



which defines the matching function Ai{y). Formally, 
I — y X is the expansion parameter so Iny = ln(l + 
(1 — y)) — > 0, which was used to simplify Eq. The 
matching condition at one-loop is therefore 



d^^ 



^ 1 



M 



(69) 



The moments of this expression are (see Appendix ^ 

,2 \ JV-i -, JV-l 



The matching condition has no large logarithms, and so 
does not depend on N in the iV — > cx) limit, except 
through the argument of a^. 



B. Breit Frame 

The computation of the matching condition in the 
Breit frame is similar to that in the target rest frame. The 
tree-level graph is given in Fig. |S1 It gives the matching 
relation Eq. H53|l where the quark distribution operator 
in the effective theory is now 



Mn{M) 



2tt 



TV — 1 TV — 1 1 poc 

_ 21n^ W 1 + 2 ^ ^ = ^J d^e— ^ [^^W^]{nzW„iz,0) 



In 



In- 



(0) 



3, 0" 7 ^2 



(75) 



which as iV — !■ cx) is 



MnIM) 



2tt 



^-21n^) (lnA^ + 7i^) 



-f (In iV + 7£;)^ + — + In^ -\ 



6 

7 ^2 



(70) 

which has the same form as Eq. (|48|l . except the external 
fields are n-collinear quarks instead of ultrasoft quarks. 
The collinear Wilson lines are required by coUinear gauge 
invariance. 

The one-loop graphs in the Breit frame are given in 
Fig. El All the graphs except Figs. I3e-g) are also 
present in the effective theory, so only these graphs 
need to be computed for the matching condition. The 
graphs are identical to the corresponding graphs in the 
target rest frame, so the matching condition is the 
^72) same, Eqs. (|68I69|) . with the replacement of Eq. H48|l by 
Eq. 
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(j) 



FIG. 9: One loop correction to the electromagnetic current product in the Breit frame. Graphs (a), (b), (e) and (h-j) also have 
mirror image graphs where the gluon attaches to the other side. 



VII. RENORMALIZATION OF PARTON 
DISTRIBUTIONS 

Parton operator renormalization gives the standard 
DGLAP evolution for parton distribution functions. The 
computation of the anomalous dimension in the target 
rest frame is identical to the computation using QCD 
quarks in Ref. jfi]. This calculation is sketched here, so 
that results can be compared with the computation in 
the Breit frame. We will only discuss renormalization of 
quark operators, which gives the evolution of the flavor 
non-singlet quark distribution. Singlet evolution mixes 
quark and gluon operators. This complication does not 
shed additional light on SCET, and will be omitted here. 



A. Target Rest Frame 

The renormalization of the (non-singlet) quark opera- 
tor can be determined by evaluating its matrix element 
in a free quark state of momentum p. The spin-averaged 
tree-level matrix element of Oq{wp'^) is given by Fig. El 
and is 



(Tp\Oq{wp'^)\p) = S{l~w). 



(76) 



The one loop graphs are those of Fig. E| as well as 
wavefunction graphs. To compute the anomalous dimen- 
sions, we only need the divergent parts of the diagrams. 



Fig. llUf a') is the discontinuity of 

g^Cp f d^fc 1_ , 1 



47r 



7 (^2^^^ 



n ■ {p — r) + 20+ 



1 



-^7 



n - {p- k-r) +iO+ '\p - fc)2 + iO+^'^'k^ + iO+ ' 

(77) 

Evaluating the k^ integral by contours, doing the kj^ 
integral, and using the substituting fc+ = zp+ gives for 
the infinite part 



h = z . „ ^ / dz 



(1-z) 



IGir^eJo {l-w + iO+){l- z-w + iO+)' 

(78) 

The integral has a discontinuity for < w < 1. As long 
as w ^ 1, the discontinuity is 

Disc/i = ^TT^- / dz{l~ z)S{l~ z~w + iO+) 



Stt^c 1 — w Jq 



in'^e 1 — w 



(0 < w < 1). 



(79) 



The singular terms as u; — *■ 1 can be evaluated by a 
method similar to that used for Eqs H64I67|I to give 



Disc Ii 



g'^Cp ( w 



'(0<w<l). (80) 
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(a) (b) (c) 

FIG. 10: One loop correction to the quark distribution function. The double line is the ultrasoft Wilson line Y . 



Figure EJb) is the discontinuity of 



9^Cf 
47r 



A'^k 1 
(2^2 



{p - kf + iO+ 

i>-^ 1 



n • (p - /c - r) + iO+ (p - ky + zO+ + iO+ ' 

(81) 



The same manipulations as the previous case give 
h = „ / dz 



so that 



Disc/2 = i 



167r-^e 



1 — z — w + jO+ ' 



(82) 



IGTT'^e Jo 



dzz(-27ri)(5 (l - z - w + iO+) 



{l-w)e{0<w < 1). 



(83) 



Figure EJc) vanishes, since rt^ = 0. Subtracting half 
the wavefunction graph for each external quark line gives 



167r2e 



5{l-w). 



(84) 



The infinite part of the matrix element is the sum of 
twice Eq. ((HOI, Eq. jHSJ, and Eq. 

g:.(o<.<i) 



w 



1 — w 



+ 1 ~ w — 2^i^ ~ ^) 



87r2e 



1 + 3,,, 

+ S{l~~w) 

(1-w) 2 



e{o<w <i) 



(85) 



in terms of the standard Altarelli-Parisi splitting kernel 



P,^g{w) = Cf 



— + -5{l - w) 



' (0 < u> < 1) . 

(86) 



Equation H85|l gives the operator renormalization equa- 
tion 



Of ( 



yp 



+\ - 



f^z(l)0A^p-), (87) 



with 



Z (z) = (5 (1 - z) + 



(z). 



(88) 



If one writes yp^ — fc+ and wp^ = then 

Of (fc+) - f^) (^+) , (89) 



which makes no reference to the quark momentum p used 
to compute the renormalization factor. The renormaliza- 
tion is invariant under boosts in the z direction, under 
which -|- components of momentum all get rescaled by a 
common factor, p+ Ap+. Eq. (|89|l is valid for fc"*" > 0. 
One can derive a similar expression for fc"*" < 0. The two 
expressions can be combined into 



Of (fc+) 



(90) 



which is valid for either sign of 

Differentiating Eq. (|90|l with respect to /i gives the 
renormalization group equation for the quark distribu- 
tion operator 



£+/k+ ^ 

(91) 



Taking moments (see Appendix ^ for the definitions) 
gives 



M [o, (fc+)] 



72,^m(,±^°°^ [O, (fc+)] , 



(92) 



where 

72, JV 



TT 

«^(a*) 

27r 



Op 



Af^ [P,^, (z)] 
^ 1 



J=2 



J iV(A^ + 1) 



1 



is the anomalous dimension. For large values of N, 
as(M), 



l2,N 



27r 



-Of [4 In TV - 3] 



(93) 



(94) 
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The infinite moments [Oq are local twist-two 
quark operators ipuft-iin ■ D) so the target 

matrix elements of [Oq (fc+)] give the familiar mo- 
ment sum rules for deep inelastic scattering. The mo- 
ments are defined by integrating over [— cxi,cxd], 
whereas the moments of quark distribution functions are 
over < a; < 1, i.e. < fc+ < P+. The matrix element 
Eq. H5U|) vanishes for |x| > 1, and its value for negative 
values of x is related to the antiquark distribution, 



fq/A-x) ^ ~fq/T{x) 



(95) 



Thus the matrix elements of M'^ for even n, which are 
the target matrix element of local twist-two operators, 
are equal to the even moments of the structure function, 
which sums over quarks and antiquarks. The matrix el- 
ements of for odd n vanish, and do not imply any 
sum rule for the odd moments of the structure function. 
The anomalous dimensions of the local twist two opera- 
tors agree with the moments of the Altarelli-Parisi kernel, 
Eq. 

In dimensional regularization, the finite parts of the 
one-loop quark distribution operator matrix element van- 
ish on-shell, since = 0, = 0, and the loop inte- 
grals are scale less. The matrix element is therefore given 
by its tree-level value 6(1 — x). The matrix elements of 
the twist-two local operators also vanish on-shell, for the 
same reason. This also shows that the quark distribu- 
tion operator is equivalent to its moments (this is not 
true for the B decay shape function 0]). The matrix ele- 
ment for off-shell quarks contains logarithms of p^//Lt^, so 
the structure function should be evolved down to some 
hadronic scale before taking the target matrix element so 
that there are no large logarithms in the matrix element. 
An important difference from the B X^j shape func- 
tion ^ |7| is that there is no additional matching that 
has to be performed at the scale Q^X^, i.e. Q'^/N'^. 

The moment analysis for deep inelastic scattering has 
been discussed in detail, even though it may be familiar 
to many readers. The reason is that there are important 
differences between the structure functions in deep in- 
elastic scattering and the shape function in the decay of 
heavy mesons, having to do with the range of the fc+ inte- 
gration in the moment of the quark distribution operator. 
The results for B decays will be presented elsewhere 



B. Breit Frame 

Structure function evolution in the Breit frame is given 
by the graphs in Fig. inia-d,h-j). Many graphs van- 
ish using = — 0. The only non-zero graphs are 
Fig.|lb,h-j). 

The ultrasoft graph. Fig. IHJb), appears to have the 
form of the ultrasoft vertex correction Eq. H27|l with an 
additional propagator for the intermediate ^„ propaga- 
tor. There is one very important difference, however. At 



scales below Q'^{l — x), the momentum of the interme- 
diate quark is of order QX, whereas the momentum of the 
ultrasoft gluon is of order QX^. Thus the power counting 
rules of the effective theory imply that the ^„ propagator 
denominator in Eq. (|27() . l/[n-{p2~k)] must be expanded 
in a power series in k, the momentum space analog of the 
multipole expansion."^ The resulting integral vanishes, so 
Fig. inib) is zero. 

The collinear graph Fig. I^i) gives 

47r 2 2^ ' n • (p - r) + iO+ 



d'^k 



1 



fin ■ (p — k) fiUa 



i-k)' 2(p - fc)2 -t- iO+ 2 k^ + iO+' 

(96) 



and Fig.intj) gives 

^ , ^ 7:Tr-p^0,7'^ 
47r 2 2^ ' ' 



n ■ (jj — r ^ k) + iO+ 



d'^fc a 1 u ?^ • (P - fc) 7^ 1 
jn r7 



{2TT)d n-k' 2ip-k)^ + i0+ 2 fc2+iO+' 

(97) 

where the relative minus sign in n • fc is because the vertex 
is from W'' rather than W. The sum of the two graphs 
has 



n ■ (p ~ r) + iO+ n ■ (p — r — k) + 20+ 
n ■ k 

[n ■ {p — r) + iO+] [n ■ {p ~ r — k) + iO+] 
and gives the integral 



(98) 



'ig'^CFT^"'p+ 
d^fc 



1 



n ■ (p — r) + iO+ 

1 n ■ (p — k) 1 

(27r)'* n ■ {p~ k ~ r) + iO+ {p- kf + iO+ k'^ + iO+ ' 

(99) 

which is the same as Ji in Eq. I|77|l , and so gives Eq. H8()|l 
for the running of the structure function. 

There is no analog of graph Fig. lTUT b). so its contribu- 
tion Eq. is missing. It is proportional to 1 — w ~ A, 
and so is of the same order as higher order terms in the 
power counting which we have dropped. 

The wavefunction contribution from the collinear 
graph is the same as the full theory wavefunction con- 
tribution, and gives Eq. H84I) . The structure function 
evolution kernel in the Breit frame is 



Pq^q{w)-CF{1 



(100) 



^ This is analogous t o th e multipole expansion for the ultrasoft 
fields in NRQCD IS ll3l . The multipole expansion is automatic 
if one makes an additional field redefinition; see the discussion 
after Eq. ITIin . 
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The moments of (1 — w) vanish as for large N. One 

can use Pq^q{w) for the evolution kernel in Breit frame, 
since the difference from Eq. 1)1 00|) is higher order in A.** 
The entire running of the structure function in the 
Breit frame is from collinear graphs. The ultrasoft graphs 
vanish. This same result has also been obtained in study- 
ing the renormalization of event shape variables using 
SCET [13. 

In SCET, one usually makes an additional field redef- 
inition 



-oo)ef (x) 



(101) 



where is an ultrasoft Wilson line. The new ri-collincar 
fields s}^^ no longer interact with ultrasoft gluons. In 
terms of these fields, Eq. H75|l becomes 



47r 



(nz) 



Yn{nz, -oo)ty„(z, 0) ^^^(0, -c5o) 



(0), 



(102) 



After the field redefinition, the collinear graphs are the 
same as before. The ultrasoft graphs have propagators 
which are independent of the collinear quark momentum 
p, since the collinear quark momentum p does not enter 
into the Wilson line. This automatically enforces the 
multipole expansion. The ultrasoft graphs vanish, since 
they are scaleless and the entire running of the structure 
function is from the collinear graphs. This is precisely 
what we found above. 



VIII. DEEP INELASTIC SCATTERING AT ONE 
LOOP 

The deep inelastic scattering structure function to one- 
loop has been computed in Ref 15]. The moments of the 
non-singlet structure function F2/{2x) are 



M 



N 



1 



-B. 



NS 



An{h), (103) 



where Ajsi^n) are the matrix elements of the twist-two 
operators renormalized at /i. The anomalous dimension 
7iv is equal to ^2,n in Eq. (|?^ . and 



^2,N — '^F 



2 ^1 ^1^1 
N{N -l-l)^j'''^^s^i 



s=l j=l 



* Consistent power counting in the Breit frame automatically im- 
plies that 1 — a;~A<^l, so the Breit frame anomalous dimension 
Eq. llOOl is not valid away from x — > 1. The target rest frame, 
however, can also be used for x far from 1. 



N N + 1 
in the MS scheme. For N 



2 

7p 



(104) 



oo, 



21n iV + 31niV- — 



(105) 



Moments of the longitudinal structure function Fl vanish 
as 1 /N as ^ cx) . 

These results can be compared with the calculations 
in this paper. The tensor structure in Eq. (|47|) shows 
that Fl = 0. The result for Fi = F2/{2x) is the product 
of the square of the matching coefficient at Q, Eq. 121|l . 
twice the running from to Q'^/N using Eq. (|37(l . the 
matching coefficient at Q'^/N, Eq. H73|l . and the running 
from Q^/n to ^ using Eq. 



M 



N 



An in) 



C'^(Q)e^^^ ln[(Q/\/]v)/Q] 

X [1 + Mn{M)] e72ln[M/(Q/V^] + o (ct2^ 



zvr 



In" 



Bin 



2 3 

2 



1 1 ^' 



0\ 



An 



'9 - TT^ /3 - 2ln N + 6\iiN 



(106) 



Using Eq. H94|l for the anomalous dimension, Eq. H106(l 
agrees with Eqs. Ill 0311 ()5t . 



IX. SCET ANOMALOUS DIMENSIONS 

There is a non-trivial relation between the anomalous 
dimensions and the matching conditions in the effective 
theory. The anomalous dimension for scaling the current 
between and Q"^ /N has the form 



7i(m) = M as{y^),\n 



(107) 



The anomalous dimension can depend on Q since that 
is a label on the SCET fields. It does not depend on TV, 
since at this stage, a single current is being renormalized, 
which has no information about 1 — x. The anomalous di- 
mension for scaling the structure function between Q"^ /N 
and n^ has the form 



72(/.) = B{as{fi)MN), 



(108) 
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where the second argument has been chosen to be IniV 
rather than N for later convenience. The anomalous di- 
mension can no longer depend on Q, since it has been 
integrated out of the effective theory. It can depend on 
N, since the parton distribution operators know about 
1 ~ X. At the matching scale between the two regimes, 
one has the matching coefficient for the structure func- 
tion 



and 



C as (a*), In 



(109) 



The matching coefficient can only depend on the relevant 
mass scale at that point, which is the invariant mass of 
the hadronic jet, Q^/iV. 

The final answer for the structure function is given by 
taking the square of the matching coefficient at Q, run- 
ning it using twice Eq. (|107|l to ~ Q^/N, multiplying 
by C (note that the matching factor is multiplicative), 
scaling by Eq. H1U8|) to a low scale, and then taking the 
target matrix element. The result must be independent 
of the choice of matching scale fi. This gives 



-M^C(^a.(^),ln^ 



72 (m) - 271 (/i). 



(110) 

Let D12C be the derivatives of C with respect to the 
first and second argument. Equation H110|l gives 

(/3(a,)i?i +2Z?2) In C = 72-271. (HI) 

where 



f3{as)C'Aas) + 2{r + l)Cr+i {as) - (116) 



if r > 1, where the prime denotes a derivative with re- 
spect to Og. 

The conditions on Ci are the usual relations that the 
logarithms in the matching condition arc the difference 
of the anomalous dimensions of the theories on either 
side. Here, however, we have the additional constraint 
2Ai = —Bi relating the two anomalous dimensions. 
Equations 137173194(1 give for the one-loop values 

Ztt 



Ai (as) 
Bo {ois) 



TT 

3a£ 
"2^ 



Cp, 



2a, 



Cf, 



Ci (as) = —^Cf, 
47r 



(117) 



and Ar = Br = 0, r > 1 and = 0, r- > 2, which satisfy 
Eq. H115|l . since P{as) starts at order a^. 



X. EXPONENTIATION 



Letting {P{as)Di + 2D2) InC = F, one gets 
2 



(112) 



F (^a„ln7V + ln^j = B (a„ In iV) - 2 A In g 

(113) 

Equation H113fl implies that A, B and F are at most lin- 
ear in their second arguments to all orders in a, . This has 
been shown previously Il6| by analyzing the momentum 
integrals of the Feynman graphs. Letting 

71 = Ai (as)ln^ + ^0 (tts) , 

72 = Bi{as)\nN + Bo{as), 



F ^ Fi(a,)ln 



+ Fo{as)., 



InC = Y,Cr{a,)W^ 



(114) 



gives 



2Ai (a,) = -Bi (a,,) = -Fi (a,) , 

Fo(a,) = Bo (as) - 2ylo (a,) , 

/3(a,)C; (a,)-f 4C2(a.) = Fi, 

P{as)C'^{as) + 2Ci{as) = Fq. (115) 



The final result for the moments of the structure func- 
tion computed using SCET is 

(118) 



with 



h - I ^ 27i(M), 



(119) 



where /ig is some reference scale of order a few GeV, 
and C{Q) M, 71 and 72 are given in Eqs. H21I74I37I94|I . 
and Ajv is the target matrix element of the local twist- 
two operator. All the large logarithms are contained in 
the exponent I — Ii + l2- Comparing with Eq. JSJ, we 
see that /o is obtained by integrating the anomalous di- 
mension using the one- loop values for Ai and Bi, fi is 
obtained by integrating the two- loop values for Ai,Bi 
and the one- loop values for Ao,Bo, etc. 

In deriving Eq. pi8|) . we have assumed that > 
Q^/N > fil. For fixed and N > Q^/^il, this inequal- 
ity no longer holds. The effective theory result for such 
large moments can be computed as follows. Match the 
current at Q and evolve using 71 from Q to fiQ. At the 
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scale ^0, one computes the matrix element of the time- 
ordered product of two currents in the target hadron. 
This is a non-perturbative computation since /xq ~ Aqcd 
is a hadronic scale. For moments with N ~ Q'^/iiq, 
the scattering process is in the resonance region since 
^ /iQ ^ ^QCD' ^^"i described by exclusive form- 
factors rather than the deep inelastic structure function. 
In equations, the result is 



(120) 



where 



— 271 (pi), 



(121) 



and A'j^{fio) is the N^^ moment of the matrix element of 
the time-ordered product of two SCET currents at the 
scale /io- 

The computation of moments is not uniformly conver- 
gent over the parameter range. One can formally take 
the limit —> 00 first, and then study arbitrarily large 
moments, since Q'^/N is now always much greater than 
Hq. However, if one takes large but finite, then for 
large enough N, Q'^/N can become smaller than /ip. The 
relevant limit for fixed target experiments is to keep the 
beam energy fixed, in which case the maximum value of 
is 2MtEx, and depends on x. 



A. Landau Pole"" 

The anomalous dimension integrals in Eq. H118|l con- 
tain as (a*) where n varies between Q and /xg. As dis- 
cussed in the preceeding paragraph, this is true even for 
N ^ 00 where formula Eq. H118|l must be modified. The 
SCET computation therefore does not suffer from any 
Landau pole singularities For large enough mo- 

ments, or equivalently, for 1 — x ~ Aqq^/Q^, the scatter- 
ing cross-section depends on resonance physics; this is a 
true non-perturbative contribution, not a breakdown of 
the perturbation series. 

The relation 2Ai = —Bi allows us to write the renor- 
malization group integration as a double integral 0. It 
is easier to derive the result starting from the answer. 
Consider the integral 



Ic = 



dy 

l/N V 



VVQ An 

2ri (a, (a^)) — + r2 {as iVyQ)) 

MO ^ 



(122) 



Changing the order of the integrals gives 



q/Vn a* 
q/Vn 



2ri {as (Ai))ln^ + 2r2 (as (a*)) 



d/i 



2ri {as (/i))lniVl 



(123) 



which has the same form as the integration of anomalous 
dimensions from Q — > Q/^/n and from Q/\/n fiQ. 
Let 



d/i 



'MO 

Comparing with Eqs. (|114I117I119|) gives 

271 = 2ri In ^ + 2r2 + Eg 
72 = 2Filn7V + r3, 



F3(a,(M)). (124) 



(125) 



so that 



Fi — —Ai — -Bi — —Cf, 

Z TT 

^3 = Bo = --^r^Cp, 
In 



(126) 



where the last entry on each line is the one-loop value. 
Eq. 1)1221) can also be written as 



/. = / dz'— 

1 — z 



dfc2 



+r2 {as {{1 - z)Q^)) 



Fl {as{{l-z)k'))^ 
(127) 



The equality between Eq. (|122|l and Eq. H127|l is derived 
in Sec. 5 of Ref. The usual expressions for the re- 
summed structure function derived using factorization 
methods have the form Eq. H122I) with fiQ /ioy^ or 
Eq. H127|l with /ig/(l — z) /lXq Q, and both suffer from 
Landau pole singularities 0| , as does the ratio of the mo- 
ments of the Drell-Yan cross-section to the square of the 
moments of the deep inelastic cross-section, A^r, which 
is independent of the hadron structure function ^ , 



■InA 



N 



- 1 



'3' dA-2 



+F2 {as {{I - z)Q^)) 



(128) 



^ I would like to thank I.Z. Rothstein for discussions on this section 



However, the corresponding SCET computation does not 
have any Landau pole singularities. The Landau pole 
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"2 




FIG. 11: Wilson line along the ni and n2 directions. 



the ultrasoft and collinear graphs are related, so that 
Ai — —Ti can be determined from the ultrasoft graph 
alone. The anomalous dimension Eq. p29|l is computed 
with two Wilson lines with nf ^ 0. The ultrasoft anoma- 
lous dimension computed earlier is obtained by taking the 
limit ni — !■ n, 712 — > n, which are both null vectors. Using 
pi (X ni and P2 oc n2, one has 



singularity arises on converting the SCET anomalous di- 
mension integration to the form Eq. (|122|l . The conver- 
sion is only valid ior / N > fj,Q in deep inelastic scatter- 
ing, and Q^/N^ > (Iq in Drell-Yan.^ The factorization 
form, Eq. (|127I128|) introduces a spurious Landau pole 
singularity, because it is equivalent to using Eq. I|118|l for 
large moments, where it is no longer valid, rather than re- 
placing it by the correct expression, Eq. (|120|l . discussed 
earlier. 

The SCET form for A^r does not have a Landau pole 
singularity, but instead contains a non-perturbative res- 
onance contribution. The resonance contributions enter 
A jv for large moments even though A jv does not depend 
on the structure function, because non-perturbative ef- 
fects enter Drell-Yan at TV ^ Q/mo before they enter deep 
inelastic scattering a.t N Q'^ / fi^, so they appear in Ajv 
beginning &i N ^ Q / ^q. 

The SCET form for the shape function in B decays 
also does not have a Landau pole singularity. 



B. Cusp Anomalous Dimension 



(2ni • 712) 



2 2 ' 

P1P2 



and 



In 77 



2 pIpI ■ 



(132) 



(133) 



The graph with p| — > in pure dimensional regulariza- 
tion replaces pf by /i^ , so that 



In = In ■ 



Comparing Eq. (|130I134I) with Eq. H114|l gives 

r, {as) = -Ai (a,) = Fi (a,) = ]-Bi (a,) 



(134) 



(135) 



Using Eq. ((nS|) for Ai, Bi, and Eq. ((TTTjl for ^0, in 
Eq. H118|l gives Fm including all terms of order as that 
do not vanish as — > 00, and including the first two 
series /o,i in Eq. 



Fi is known as the cusp anomalous dimension, because 
it can be computed from the anomalous dimension of a 
Wilson line with a cusp m (seeFig.inj- The anomalous 
dimension of the Wilson line in Fig. ^] has the form 



(2wi ■ 712) 

2 2 



(129) 



since it must be Lorentz invariant, and invariant under 
rescaling of the vectors 77^. As 77 —> 00, 



iw {oLs,vi) In 77 Tc (as) , 



(130) 



which defines the cusp anomalous dimension Fc p^ . To 
two loop order T6j 



a. 



Cf + 



V TT 



Ca 



67 
36 



TT 

12 



-TpNp 



(131) 



XI. CONCLUSIONS 

The deep inelastic structure function was computed in 
the a; — > 1 endpoint region using SCET. The effective 
theory properly separates the scale Q^, Q^{1 — x) and 
-^QCD' ^^'^ allows the problem to be analyzed one scale 
at a time. The structure function is equivalent to local 
twist- two operators below the scale Q'^{1 — x). The scale 
Q^{1 — x)"^ does not play a special role in the analysis 
of structure functions. This is very different from the 
case oi B ^ Xgj decay, where the shape function has an 
anomalous dimension between Q^{1 — x) and Q^{1 — x)^, 
and stops running below Q^{1 — x)^ 130 ■ Consistency 
of the effective theory implies that the SCET anomalous 
dimensions are linear in In /i to all orders in perturbation 
theory. The SCET formulation also avoids the Landau 
pole singularity in the resummcd cross-section. 



We have seen earlier that the equality of the anomalous 
dimensions in the target rest frame and Brcit frame im- 
plies that the In fi terms in the anomalous dimension from 



The extra A'^ for Drell-Yan arises because of the difference in 
kinematics. 
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APPENDIX A: MOMENTS 

Let f{z), g{z) and h{z) be functions defined for < 
z < 1. The convolution of / and g,h~f-kg~g-kf is 
defined by 

h{z) = / da; / Ay 5{z - xy)J [x) g {y) 



The moments of a function are defined by 



MN[f{z)] - 



and satisfy the relation 

MNif^g] = [/] Afjv [3] . 
The moments we need are 



Mn[S{1~z)] = 1, 
1 

TV' 



Mn [1] 
Mn [1 - z] 



1 



1 - z 



N{N + 1)' 

N-l+r ^ 

= Hr — HN^i+r, 



(Al) 



(A2) 



(A3) 



M 



N 



/ln(l - z) 



N-l 



V 1 - ^ 

where the harmonic number Hr is defined by 

1 



(A4) 



(A5) 



The large N limits of the moments are (N = Ne'"^) 



Mn [5{l - z)] 
Mn [1] 
Mn [1 - z] 



Mn 



M 



N 



1 - z 

ln(l - z) 

1- z 



1, 
0, 
0, 

Hr - In TV, 



lln^N+lc{2). (A6) 



For functions defined on (— cx),oo), one can define the 
moments over [0, 1] as in Eq. ljA2p as well as the half- 
infinite and infinite moments 



M+[f{z)] = / dzz^-V(^), 



M^[f{z)] = / dzz^-VW, 

/OO 
dz z^-'f{z) 
-00 

= Mj^[/(z)]+M+[/(z)]. (A7) 
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